Stable two-dimensional ferromagnets made of regular single-layered lattices of 
single-molecule nanomagnets on substrates 
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We propose that stable two-dimensional (2D) ferromagnets can be made of regular single-layered 
lattices of single-molecule nanomagnets with enough uniaxial magnetic anisotropy on appropriate 
substrates by controlling the inter-nanomagnet magnetic interaction. Our Monte Carlo simulated 
results show that such ideal 2D ferromagnets are thermodynamically stable when the anisotropy is 
strong enough. If the anisotropy energy equals 80 K, approximately that of the Mni2, the Tc varies 
from zero to 15 K depending on different inter-nanomagnet coupling constants. Such stable spin 
systems, experimentally accessible, should be promising for information applications. 

PACS numbers: 75.10.-b, 75.75.+a, 75.30.-m, 75.90.-|-w, 81.05.Zx 



Single-molecule (SM) nanomagnets attract huge in- 
terest because they show quantum tunnelling, interfer- 
ence, and coherence effects at ultra-low temperatures 
and have various potential applications, such as mag- 
netic information storage and quantum computation[l), 
i, i, i, [i, @, 0, H]. An SM nanomagnet usu- 
ally has large spin 5* and strong uniaxial anisotropy, 
and there is an easy axis for spin orientation[8]. 
Hence it must overcome an effective energy barrier 
Ucs to achieve a spin reversal at low enough tempera- 
ture. For example, the famous Mni2 SM nanomagnet 
([Mni20i2(02CCH2Br)i6(H20)4]-4CH2Cl2) has S = 10 
and its spin reversal energy barrier Ucff is upto 74 k[8|,[^. 
The record-breaking value of Ucs is 86.4 K in the case of 
[Mnf^02(Et-sao)6(02CPh(Me)2)2-(EtOH)6], which has 
S = 12[3- 

SM magnets can be made to form single isolated mag- 
nets, dimers, and other clusters[l|, [3, H, 0, HI- Through 
chemical methods, they can be synthesized in the form 
of three-dimensional (3D) crystalline structures and two- 
dimensional (2D) networksji, E [Hi [H, [11]. The ba- 
sic building blocks in these structures are single mag- 
nets. The effective inter-magnet magnetic interactions 
usually are antiferromagnetic (AF) and cannot be eas- 
ily controlled. Thus, it is quite difficult to make stable 
ferromagnets from such systems of SM magnets. Fortu- 
nately, single layers of SM nanomagnets have been re- 
cently achieved on gold and silicon surfaces 14,15, [l^ . 
This breakthrough makes us believe that stable ideal 2D 
ferromagnets can be made from regular 2D lattices of SM 
nanomagnets, although Mermin- Wagner theorem allows 
no finite-temperature phase transition in 2D isotropic 
Heisenberg spin systems 171, Il8|. 

Here we show that stable 2D ferromagnets can be 
made of regular single-layered lattices of SM nanomag- 
nets with enough uniaxial anisotropy on appropriate sub- 
strates by controlling the inter-nanomagnet magnetic in- 
teraction. The spin of the nanomagnet is described as a 
large spin. The easy axis is perpendicular to the 2D spin 



lattice. We describe the 2D ferromagnets by a 2D quan- 
tum Heisenberg model with the uniaxial anisotropy and 
relatively weak ferromagnetic (FM) inter-spin coupling. 
This model is reasonable and liable because the dipo- 
lar spin interaction almost is reduced to a purely AF one 
and the inter-spin coupling can be controlled by choosing 
right substrates and manipulating the inter-nanomagnet 
distances [l9j|. We reduce the quantum spin operators to 
classical variables because both the spin value and the 
anisotropy energy are large enough [sl [ol [lo| . and then use 
Monte Carlo method to simulate the many-body physi- 
cal properties with various parameters j2Q, i21||. Our re- 
sults prove that stable 2D ferromagnets can be made of 
such SM nanomagnets for experimentally-accessible pa- 
rameters. More detailed results will be presented in the 
following. 

SM nanomagnets such as the Mni2 have large spins and 
their inter-spin exchange interactions usually are very 

weak compared to their anisotropy energies d, u, [3, [U- 

When put on appropriate surfaces, such SM nanomag- 
nets can have the same uniaxial anisotropy energy and 
adjustable inter-spin interactions [l3. [isl \w^. Thus, the 
2D ferromagnets can be described by the Hamiltonian 



i/ = -K^(5f)2-^J,,S',.5, 



(1) 



where K describes the uniaxial anisotropy with the easy 
axis {z direction) perpendicular to the 2D lattice, the 
Si is the spin operator at site i in the lattice, and Jij 
is the inter-spin coupling constant between Si and Sj. 
This perpendicular anisotropy and the 2D lattice struc- 
ture make the dipolar spin interaction almost reduced to 
an AF one, and thus Jij can varies from AF to FM in- 
teraction. Here we are interested mainly in the FM case. 
The latter summation is over spin pairs between i and j. 
Because the spin 5* is very large and we are interested in 
the equilibrium properties, we can reasonably treat the 
quantum spin operator Si as a classical vector Ssi, where 
Si is a unit vector. 
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We investigate the equilibrium properties of this model 

Both Metropo- 
2i. We take an 



using Monte Carlo (MC) method[2C 




22 



lis and cluster algorithms are used 
Lx L square lattice and limit Jij to nonzero J (> 0) only 
for the nearest spin pairs without losing main physics. As 
usual, we use a periodic boundary condition. For conve- 
nience, we use a special unit that the Boltzmann constant 
fcs is set to 1. We use the parameter D = K/J to de- 
scribe the strength of the anisotropy energy relative to 
the exchange constant. As we prove, the system is al- 
ready in equihbrium after 100,000 MC steps (MCS). The 
average magnetization M, the fourth-order Binder's cu- 
mulant /74, the specific heat C, and the magnetic suscep- 
tibility X a'l'e calculated by averaging over 50,000 MCS af- 
ter thermodynamical equilibrium is reached [iSIIil. Ten 
copies of such average values are used to calculate the 
final results. 
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FIG. 1: (color online). Normalized magnetizations M as func- 
tions of temperature T for given anisotropy parameters from 
Z)=0.05 to 500. 



In Figure 1 we present average magnetizations M as 
functions of temperature T for 16 D values from 0.05 
to 500. We have M = 1 at T = because the mag- 
netization M is normalized to 1. The magnetization M 
for all the higher temperatures increases with D. There 
is a spinodal point in the specific magnetization curve 
for a given D value. When D becomes larger, the spin- 
odal point shifts rightward to higher temperature, which 
means that the phase-transition temperature increases 
with increasing D. In order to determine more accu- 
rately, we actually use the unique intersection point of the 
various U4 — T curves of different L values. This method 



is the best approach to obtain accurate phase-transition 
temperature 20|, 



21|. 



The calculated Tc as function of 
D is presented in Figure 2. In order to emphasize the 
Tc of large D values, we use logarithmic scale in the 
main plot and linear scale in the inset. Our results show 
that Tc increases with D, converging to 2.269, the Ising 



limit 2J, |25|, when D is infinite. For D = we obtain 
Tc = 0, which is consistent with established conclusion 
that there is no finite-temperature phase-transition for 
2D isotropic model 17 , l^l • A nonzero D is necessary to 
a finite Tc- 
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FIG. 2: (color online). The D dependence of Curie temper- 
ature Tc with D from 0.05 to 500. The dash line indicates 
Tc=2.269, the Ising limit. The inset shows Tc as a function 
of D in linear scale. 



We describe the deviation of the normalized spin vector 
from the easy axis by the parameter d,„ defined as 



d™-l-((s^f)/((.s^)), 



(2) 



where s = {1/L'^) J2 and (^4) is the average value of A. 
We present dm{T) for various D values in Figure 3. It is 
clear that dm is zero or tiny at low temperature. When 
D approaches to zero, we obtain d,„ < 2/3. For a given 
Z?, dm increases with T, reaching a maximum dmc at Tc. 
dmc decreases with increasing D, being already less than 
0.02 for D = 30. For D > 100, dm is actually zero in the 
temperature window in Figure 3. Indeed, when D is very 
large, the Hamiltonian is similar to the 2D Ising model, 
but it is equivalent to the 2D Ising model only when D 
is infinite. 

It is established that the inter-spin interaction can 
be changed experimentally and even adjusted intention- 
ally for special purposes [l9j. This means that J can be 
changed actually from a few K to zero. Considering that 
the parameter K in Eq. (1) can be nearly 80 K 10|, we 
conclude that the anisotropy parameter D can be 10 or 
much larger. Accordingly, Tc is 1.8 ~ 2. 269 J. With 
if = 80 K, Tc equals 2.2 K if J 1 K, and 14.4 K if 
J ~ 8 K. Because the 2D ferromagnets are in the regime 
of large D, the spin deviation dm is very small, smaller 
than 0.001, when T is at Tc, and can be considered zero 
when T is substantially below Tc- Therefore, the spin is 
nearly in the Ising limit in well-ordered FM states and 
such spin systems should be stable enough to carry in- 
formation. 
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FIG. 3: (color online). The temperature dependence of the 
deviation dm of the spin vector from the easy axis for various 
D values (from 0.05 to 30.0). 



In summary, we propose that stable 2D ferromag- 
nets can be made of regular single-layered lattices 
of SM nanomagnets with enough uniaxial magnetic 
anisotropy on appropriate substrates by controlling the 
inter-nanomagnet magnetic interaction. The quantum 
spin operator for them can be naturally reduced to a 
classical spin variable because the spin value and the 
anisotropy energy are so large. Our Monte Carlo sim- 
ulated results show that such ideal 2D ferromagnets are 
stable because the deviation of the spins from the easy 
axis is smaller than 0.1% in the ferromagnetic state when 
the anisotropy is strong enough. If the anisotropy energy 
equals 80 K, approximately that of the famous Mni2, 
the Tc varies from zero to ^--^15 K according to the dif- 
ferent inter-spin coupling constants. Our 2D model is 
reasonable and reliable because the dipolar inter-spin in- 
teraction almost can be reduced to an AF one due to the 
perpendicular easy axis and the 2D spin lattice. There- 
fore, such stable spin systems should be promising for 
information applications. 
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